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Muons produced in Extensive Air Showers (EAS) generate ring-like images in Imaging Atmo¬ 
spheric Cherenkov Telescopes when travelling near parallel to the optical axis. Erom geometrical 
parameters of these images, the absolute amount of light emitted may be calculated analytically. 
Comparing the amount of light recorded in these images to expectation is a well established tech¬ 
nique for telescope optical efficiency calibration. However, this calculation is usually performed 
under the assumption of an approximately circular telescope mirror. The H.E.S.S. experiment 
entered its second phase in 2012, with the addition of a fifth telescope with a non-circular 600 m^ 
mirror. Due to the differing mirror shape of this telescope to the original four H.E.S.S. telescopes, 
adaptations to the standard muon calibration were required. We present a generalised muon cal¬ 
ibration procedure, adaptable to telescopes of differing shapes and sizes, and demonstrate its 
performance on the H.E.S.S. II array. 
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1. Introduction 

Charged particles within Extensive Air Showers (EAS), generated by energetic /-rays or Cosmic 
Rays penetrating the Earth’s atmosphere, travel faster than the local speed of light, thereby produc¬ 
ing Cherenkov radiation, subsequently detectable by Imaging Atmospheric Cherenkov Telescopes 
(lACTs). In order to record snapshots of EAS, fast imaging cameras are employed, with typical in¬ 
tegration times of 16 ns. The nature and energy of the primary particle may be determined from the 
shape and intensity in photoelectrons of the recorded images respectively. The combined optical 
efficiency of the various telescope and camera components directly influences the image intensity 
and can significantly affect the energy reconstruction. 

Muons produced in hadronic EAS are more penetrating than lighter leptons, and may scatter to 
travel further from the shower axis; this reduces the amount of radiation from parent hadronic EAS 
recorded in the same image as a muon. The ring-shaped images generated in the camera are easily 
identifiable; the use of these images for optical efficiency calibration is a well-established method 
within the lACT community [|I]]. Comparing the recorded image amplitude for a given muon ring 
to the theoretical prediction provides a measure of the light collection efficiency of the telescope 
system. The mean of this parameter distribution over a given time period is termed the muon 
efficiency, e^. 

Within the analysis chain, this muon efficiency is used to correct energy estimates by a factor 
accounting for the discrepancy between the current state of the system and the relevant Monte 
Carlo (MC) reference state. The optical efficiency can also have a large effect on the Instrument 
Response Eunctions which must be accounted for by regular calibration. 

H.E.S.S. (High Energy Stereoscopic System) is an array of lACTs situated at 1800m altitude in the 
Khomas Highlands in Namibia |^], enhanced in 2012 by the addition of a fifth telescope of 600m^ 
mirror area (hereafter CT5). This enables the energy range of H.E.S.S. to be extended at the low 
end, pushing down the energy threshold to ~ O(lOGeV). Whilst the muon calibration technique 
has been utilised for telescopes of the H.E.S.S. array since its inception in 2003 [||, standard 
assumptions of the muon reconstruction procedure, in particular a circular mirror shape, were no 
longer valid for CT5. Here, the adaptations subsequently made to the standard muon calibration 
are outlined, which render the algorithm suitable to varying lACT mirror configurations. Improved 
performance over previous implementations of this calibration procedure is also demonstrated. 

2. Muon Calibration Procedure 
2.1 Muon Image Identification 

The raw image is initially cleaned using tail-cuts (dual-threshold: all pixels above an upper thresh¬ 
old and neighbouring pixels above a lower threshold are kept), with upper and lower thr'esholds of 
7 and 5 photoelectrons respectively. An image is subsequently classified as ring-like if a circle is 
found using the Chaudhuri-Kundu formula [Q]. In this case, the circle fit is re-performed using the 
results of the first fit as seed parameters and excluding pixels lying further than 1.5 pixel widths 
away from the seed ring radius. This helps to refine the ring parameters (radius, centre position, 
ring width) by reducing fit biases due to scattered light from parent hadronic showers. Intensity 
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Figure 1: Line integration along the mirror surface as a function of azimuth angle for a given impact position, 
shown for CT5 with the true mirror profile (solid line) and under a circular approximation (dashed line). 


information for “on ring” pixels (pixels within 0.5° radial distance either side of the fitted ring 
radius) is then taken from an uncleaned version of the image. 

A set of cuts optimised for this algorithm and for different telescope types, are used to select muon 
images for calibration. These cuts ensure that the selected ring images are (i) well above threshold, 
(ii) well contained, (hi) ring-like and (iv) that the influence of hardware issues is minimal. Pixel 
related cuts are made on the total number, Npix (reason i), the number in the outermost layer at the 
camera edge, NedgePix (reason ii), the average number of neighbouring pixels, (NN) (reason hi), 
and the number classified as broken, (reason iv). Additional cuts are made on the ring radius, 
r (reason i) and on the outermost radius of the ring, r-\- Pr (reason ii), where Pr is the ring impact 
parameter (the distance from the centre of the ring to the centre of the camera). 

2.2 Predicted Image Intensity 

The total intensity Ipe of Cherenkov light produced by a muon as a function of azimuthal angle 0 
and impact parameter p, the distance from the muon impact position on the miiTor to the centre of 
the mirror, is given by: 


Ipe = £gI{0c,p,^,O)) = ^1 ^m{2Bc)D{p ( 2 . 1 ) 

where ft) is the angular camera pixel size; a the fine structure constant; and the muon efficiency, 
a global scaling parameter relating the total light intensity to the expectation, determined by a 
log-likelihood based minimisation procedure. The Cherenkov opening angle of light emission, dc, 
corresponds to the ring radius, r as measured in degrees in the camera frame; and D{p,(p) is the 
distance across the miiTor along which the collected light must be integrated. The intrinsic intensity 
/ = / A^^dA is integrated over optical wavelengths 300nm < A < 600nm, where the atmospheric 
refractive index is assumed to be A independent ||^. 
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Previously, D{p,^) was determined analytically assuming a circular min~or dish for H.E.S.S. I; an 
approximation no longer adequate for describing the CT5 miiTor (figure |I|). Instead, D{p,^) is 
obtained by interpolating the distance between a set of points describing the outside of the dish. 
Figure |I] shows that the variation in distance with azimuth angle is more accurately described by 
the interpolation, avoiding biases due to a circular approximation, such as an apparent shift in 
azimuthal position of the central hole. The azimuthal muon ring intensity profile ( ^T| ) lakes Ihe 
same shape as lhaf shown in figure [I] for D{p,(p) since all olher paramefers are (j) independenf. 


2.3 Muon Efficiency Determination with a 2D Fit 


The raw pixel-wise intensities for all “on ring” pixels are summed along each azimuthal direction 
around the ring, producing a ID distribution of intensity as a function of azimuth angle, I{(p). This 
raw intensity distribution is smoothed around the ring with a four pixel-widths moving average 
sliding window. After smoothing, the ID intensity distribution is broadened radially, according to 
a Gaussian distribution centered on the ring location, and spread by the fitted ring width. In all 
stages, the total integral chai'ge is preserved. The 2D smoothed image is fitted with a 2D model 
using a pixelwise log-likelihood, with the overall efficiency, impact position on the mirror and ring 
width as free parameters of the fit (the ring centre and ring radius are kept as fixed variables). In 
contrast, the algorithm previously implemented used a ID fit to the azimuthal ring intensity 
profile to determine e^. The pixel likelihood is calculated from a Gaussian approximation after |§], 
for the case of high model prediction compared to the observed signal. An example muon event 
recorded by CT5 and the 2D model expectation is shown in figure §. 

Events passing further cuts on the impact parameter (distance from the mirror centre) and on the 


ring width, (for reasons ii and iii introduced in section |2.1[ ) are retained for the calibration. The 
telescope-wise muon efficiency is given by the mean of the distribution from single events per run 
(approx, half an hour of observations). The obtained muon efficiency is compared to a reference 



Figure 2: Muon event as seen in the CT5 camera (left) and the expectation from the 2D fit (right). Colour 
scale corresponds to the charge in photoelectrons per pixel after calibration. Statistical fluctuations may 
yield negative values after baseline subtraction. 
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Figure 3: Linearity of previous and new algorithms for both the H.E.S.S. I telescopes (left) and for CT5 
(right). In both cases, an improvement over the performance of the previous algorithm is seen. 


found by performing the muon analysis on MC. This reference MC uses input telescope optical 
efficiencies corresponding to the set of lookup tables used for the energy reconstruction. The ratio 
of the reconstructed MC muon efficiency fo fhe muon efficiency obfained on dafa defines an optical 
efficiency correction facfor, applied fo energy esfimates from individual felescopes af analysis level. 

3. Performance 

Two major changes fo fhe muon calibrafion algorifhm in fhe H.E.S.S. analysis chain, fhe arbifrary 
mirror description and fhe 2D pixel log-likelihood tiffing, have led fo demonsfrable improvemenfs 
in fhe performance. Ideally, fhe correcfion facfor should be minimally biased by efficiency degrada- 
fion. This can be checked by performing fhe muon calibrafion on MC of various (known) felescope 
degradafions. The correcfion facfor found from performing fhe muon calibrafion on MC, e™"’/ 
should be linear wifh fhe rafio of fhe known inpuf felescope optical degradations. Here, e™™ is fhe 
muon efficiency obfained from using fhe calibrafion algorifhm on MC of fhe felescopes af fheir 
nominal efficiency (defined as 100%), whilsf is fhe muon efficiency obfained from using fhe 
calibrafion algorifhm on MC of fhe felescopes af a degraded efficiency. The improvemenf in linear 
behaviour wifh fhe updated algorithm is clearly demonstrated in figure |3[ 

This generic algorifhm is currenfly used on all five H.E.S.S. felescopes as par! of fhe sfandard 
calibrafion chain. The evolution of fhe muon efficiency of CT5 and one of fhe H.E.S.S. I felescopes 
during 2014 is shown in arbifrary unifs in figure Alfhough fhe overall efficiency evolufion is 
mosfly sfable wifh a gradual degradation, fhe runwise efficiencies are still subjecf fo some sfalislical 
variation. The merging of muon efficiencies over longer lime periods, such as nighlly, may reduce 
sfalislical varialion fo provide more sfable eslimafes. Disconlinuous changes are discernible in bolh 
telescopes af run number 95000 due fo camera hardware adjuslmenl, and in CT3 af run number 
100700, due fo cleaning of some camera componenls. 
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Figure 4: Evolution of the muon efficiency (in arbitrary units) of CT5 and one of the H.E.S.S. I telescopes 
during 2014. Discontinuous changes due to camera hardware adjustment in both telescopes, and cleaning of 
CT3 camera components are marked by dashed and dotted lines respectively. 


4. Conclusions 

A generic muon algorithm has been implemented for H.E.S.S. II, with the flat distance across a 
telescope mirror determined by interpolation between a set of points defining the mirror edges. 
Smooth (circular or elliptical) mirror shapes are definable with a large number of points (unre¬ 
stricted in principle). A 2D pixel log-likelihood fitting procedure replaces the previous ID 'x} fit to 
the azimuthal muon intensity profile. These alterations have led to demonstrable improvements in 
the linearity of the correction factor with optical degradation. 

The development of flexible algorithms for telescopes of varying specifications is of increasing 
importance with the advent of the Cherenkov Telescope An'ay (CTA) [^. We foresee no com¬ 
plications arising from the adaptation of this algorithm for arbitrary miiTor shapes, whilst minor 
adjustments may be required for telescopes with Schwarzschild-Couder (dual min'or) optical de¬ 
sign. 
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